Surface modes and vortex formation in dilute Bose-Einstein condensates 

at finite temperatures 
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The surface mode spectrum is computed self-consistently for dilute Bose-Einstein condensates, 
providing the temperature dependence of the surface mode induced vortex nucleation frequency. 
Both the thermodynamic critical frequency for vortex stability and the nucleation frequency implied 
by the surface excitations increase as the critical condensation temperature is approached from be- 
low. The multipolarity of the destabilizing surface excitation decreases with increasing temperature. 
The computed finite-temperature critical frequencies support the experimental observations and the 
zero- temperature calculations for vortex nucleation. 
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Since it became feasible to create and observe singly 
quantized vortices D and vortex lattices Q in di- 
lute Bose-Einstein condensates, a desire for quantitative 
understanding how vortices penetrate these unique quan- 
tum systems emerged. When the system is externally 
rotated above a certain critical nucleation frequency 
vortices enter the condensate from its boundary. How- 
ever, the exact dependence of Vli, on the external vari- 
ables and the experimental conditions remain the topic 
of continuing discussions. The reported minimum ratios 
between the external rotation and the harmonic trap- 
ping frequencies for which vortices have been observed 
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thus far include approximately 0.4 
oblate, and 0.1 §, 0.3 §, and 0.7 g _ 
geometries. Also, a number of theoretical interpreta- 
tions have been proposed for explaining these frequencies 
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First, there exists an angular frequency above which 
the free energy of an axially symmetric vortex state be- 
comes lower than that of the vortex-free state. In most 
experiments, this thermodynamic critical frequency Qc, 
has been found to be substantially below the lowest fre- 
quency for which vortices have been observed. This may 
be explained in terms of an entry barrier at the surface 
of the condensate, preventing the transition to the lower 
free energy state. This, in turn, implies hysteresis effects 
when the rotating drive is ramped above the nucleation 
threshold and back again. Hence, vortices should occur 
in condensates at much lower rotation frequencies than 
those used to create them. In fact, according to the finite- 
temperature Popov, Gl and G2 approximations, axisym- 
metric singly quantized vortex states are predicted to be 
locally energetically (meta)stable even below Qc ||2^, 
(see, however, the discussion in Ref. p6| ). 

Secondly, there exist nodeless surface excitations, or 
surfons, in the spectrum of the Bose-Einstein conden- 
sate that are localized near the surface of the condensate. 
These 'quantum tidal waves' pTf are often interpreted as 
being responsible for the instability of the condensate, 
leading to the onset of vortex nucleation as they rep- 
resent a density perturbation on the condensate surface 



that may evolve into a quantized vortex line, thus chang- 
ing the topology of the condensate. For rotationally in- 
variant systems, an analog of the Landau criterion 
may be expressed for the surfons as 



lis = mm < — 
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where huji is the excitation energy of a surfon carrying 
the angular momentum HI. This angular frequency of the 
rotating perturbation corresponds to the instability of 
the vortex-free condensate since beyond it there emerges 
anomalous negative-energy excitations with high multi- 
polarities in the spectrum of the condensate ||T2l. How- 
ever, the value for the nucleation frequency predicted by 
Eq. (0) is still too low for serving to explain the ob- 
served vortex nucleation thresholds in Refs. ^, 
Accordingly, it has been suggested that in those cases 
the external perturbation used to rotate the system has 
mainly driven the (/ = 2) quadrupole surface mode, 
for which uJi/l yields a much higher critical angular fre- 
quency. Also, the bending of vortex lines and anomalous 
vortex core modes have been argued to explain certain 
aspects of the vortex nucleation process [T^ . 

In this paper, we study finite-temperature effects on 
vortex nucleation and the thermodynamic stability of 
vortices in dilute Bose-Einstein condensates. The tem- 
perature dependence of the surfon-induced vortex nucle- 
ation frequencies, as well as the thermodynamic criti- 
cal frequency, are computed. The numerical results ob- 
tained are consistent with the experimental observations 
and support the present theoretical understanding of the 
vortex nucleation process based on zero-temperature field 
theories. Although our results are in qualitative agree- 
ment with those of Ref. |l5j, there exist quantitative 
differences besides the system parameters, since we do 
not discard any low-energy quasiparticle excitations from 
the self-consistent sums, and we do employ the local 
density approximation taking properly into account also 
the high-lying quasiparticle states. Furthermore, due to 
the thermal gas, there is a slight difference between the 
prediction of Eq. (pi) and the frequency at which there 
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emerges anomalous surface mode(s) in the spectrum cal- 
culated in the rotating frame of reference. 

The stationary generalized Gross-Pitaevskii equation 
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Ho + .g|0(r)P + 2.gn(r) -r!L, 
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for the macroscopic wavefunction (/)(r) models the Bose- 
Einstcin condensate at finite temperatures. Above, fi 
stands for the chemical potential, Ho = — (fi.^/2M)V^ + 
V{r), and the interaction strength g = Airh^a/M ^ where 
a is the scattering length for atoms of mass M. The wave- 
function is normalized according to /[|0(r)p -|-n(r)]dr — 
N and the system is rotated at the angular velocity H, 
providing the N trapped atoms with the angular momen- 
tum L. Furthermore, the density of the noncondensate 
atoms ri(r) is determined self-consistently via the relation 
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:(K(r)p-t-|^,(r)pK + K(r)|^ 
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where Uq — (e 



— (^^Eg/kBT _ is the Bose-Einstein 

distribution function and the quasiparticle amplitudes 
Uq(r),Vq{r), and eigenenergies Eq are obtained from the 
Bogoliubov equations 

Cuq{r)+g(l)'^{r)vq{r) = EqUqiv), (4a) 

C*Vq{v)+g<j,*\Y)Uq{v) = -EqVq{v). (4b) 

Within the Popov approximation employed here, C = 
7^0— M+25[|0(r)p+7T,(r)]— ili^, and q labels the quasipar- 
ticle states whose amplitudes must obey the normaliza- 
tion condition /[|Mg(r)p — |w^'(r)p]dr = Sqqi. Negative- 
energy quasiparticle excitations with positive norm are 
referred to as anomalous modes. The free energy of the 
system, F = (Hcs) — TS, is evaluated using the equation 



9 



\(j){v)\'^Av~2g / n(r)|</.(r)pdr 



+ / Eq\vq[v)\^Av (5) 

- k^T'y^\{l + nq) ln(l + Ug) - Ug Inn^]. 



We shall consider cylindrically symmetric states, for 
which (/)(r) — ^(r)e""^, the integer winding number m 
determining the number of circulation quanta in the vor- 
tex. Consequently, the equations may be separated into 
the radial, angular and axial parts, and the surfons are 
uniquely labeled by their angular momentum quantum 
number I (the radial and axial quantum numbers for 
these states equal zero). The condensate is radially con- 
fined by a harmonic trapping potential V{r) — ^Muj^r^. 
In the axial direction we impose periodic boundary con- 
ditions. The relevant physics of this system is expected 
to closest resemble spherical/oblate condensate geome- 
tries, as the comparison of the critical frequencies with 



those in Ref. [|2| also sug gest s. In order to guarantee 
the accuracy of our results |30|| , we compute the discrete 
quasiparticle states using Eqs. (^ up to Eg = SOftw^ and 
employ the local density approximation ^ for the 
remaining states. In the results presented, ujr = x 10 
Hz and the effective gas parameter Na/aho ~ 90, where 
flho = y^h/mujr denotes the harmonic-oscillator length. 

Energies per angular momentum of the surface modes 
with different multipolarities are plotted in Fig. |l| for 
a range of temperatures. The solid line is the hydro- 
dynamic prediction for the surfon dispersion rela- 
tion [Q. At finite temperatures, Kohn's theorem for 
parabolic confinement, according to which the dipole 
mode (/ = 1) should oscillate exactly at the trap fre- 
quency (jJr, is not satisfied due to the neglect of the non- 
condensate dynamics. For all temperatures, the nucle- 
ation frequency cui/l implied by the quadrupole mode is 
found to be close to the value 0.7 observed in the experi- 
ments for which the quadrupolar drive may be expected 
to be the primary source for nucleation p, ||, ^ . This 
confirms the negative result found in Ref. B for the sig- 
nificant temperature dependence of the vortex nucleation 
induced by a quadrupolar drive. That particular nucle- 
ation threshold is also expected to be fairly insensitive to 
the particle number and the geometry of the trap . 

However, the energies of the surface modes with higher 
multipolarities progressively increase at higher tempera- 
tures, and the lowest critical frequency Qs due to the 
surfon excitation is found to increase approximately from 
0.4 ujr to 0.5 LOr with increasing temperature. Also, the 
corresponding surface mode shifts to a lower value of an- 
gular momentum, see the inset in Fig. |l|, which is a direct 
consequence of the shrinking of the radius of the conden- 
sate. 

The computed radial distributions for the densities of 
the condensate, noncondensate and the destabilizing sur- 
face mode are depicted for 51 = in Fig. ^ for two differ- 
ent temperatures. These surfons are localized near the 
edge of the condensate, of which only a part is seen within 
the area of the figure. The shift in the multipolarity of 
the destabilizing surfon as a function of temperature is il- 
lustrated in the figure. As the temperature increases the 
radial extent of the condensate cloud shrinks. Hence, the 
angular momentum of the surfon localized in the vicinity 
of the condensate boundary has the larger angular mo- 
mentum the lower the temperature is. Notice the non- 
condensate bulge, due to the mutual repulsion between 
the condensate and noncondensate particles, on the outer 
edge of the condensate. 

In Fig. ||we have plotted the lowest surfon-related (o) 
nucleation frequencies Qg and two different thermody- 
namic critical frequencies. Also shown is the critical con- 
densation temperature T^ifl) = T^{1 - fl^/cj^Y^^ || 
and the nucleation frequency suggested by the anomalous 

vortex core mode ( ), calculated within the nonselfcon- 

sistent Bogoliubov approximation. Both the thermody- 
namic and the surfon critical frequencies remain nearly 
constant below Tc/2. At higher temperatures, espe- 
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cially the thermodynamic equihbrium critical frequency 
ric (•) strongly deviates from its zero-temperature value, 
in agreement with the calculation of Ref. . Further- 
more, the critical frequency, given by the energy differ- 
ence between the vortex and nonvortex states in the ab- 
sence of rotation {il = 0), is shown (■). Physically, this 
corresponds to the situation in which the noncondensate 
stays at rest, which has probably been the case in most 
of the available experiments, excluding those of Ref. 
The decrease of the critical frequency with increasing 
temperature is probably due to the neglect of the non- 
condensate dynamics by the Popov approximation whose 
predictions also for the excitation energies are known to 
deviate from the experimental values for temperatures 
T > Tc/2. 

We have also computed the surfon spectrum for an ax- 
isymmetric vortex state and found that the correspond- 
ing nucleation frequencies are only slightly higher than 
those presented in Fig. This is in accordance with 
the fact that the angular rotation frequency needed for 
nucleating vortices grows in proportion to the number 
of vortices present in the system. Moreover, a presence 
of multiply quantized vortices in the harmonic trap | |3^ 
should not affect the qualitative process of vortex forma- 
tion at the condensate boundary. 

In conclusion, we have computed the surface mode 
spectrum and the thermodynamic phase diagram for vor- 
tex stability in the T — plane for dilute Bose-Einstein 



condensates at finite temperatures. According to the 
current understanding, the highest observed nucleation 
thresholds with fl/tUr ~ 0.7 can be explained in terms of 
the excitation of the quadrupole surface modes, whereas 
the lower values 0.3— 0.4 lie close to the minima oiuji/l for 
the corresponding systems. This scenario for vortex nu- 
cleation in dilute Bose-Einstein condensates is suggested 
to be valid by our computations also at finite temper- 
atures. However, the nucleation value 0.1 coinciding 
with the thermodynamic critical frequency, cannot solely 
be explained by the present analysis. A remedy of this 
observation in the favor of the surfon instability mech- 
anism could be obtained by noting that the local fluid 
velocity around the small stirrer used in the experiment 
may well exceed the velocity of the stirrer and thus the 
effective rotation frequency could be locally much higher 
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FIG. 1: Temperature dependence of the surface mode energies 
divided by the corresponding angular momenta as functions 
of the angular momentum quantum number /. The solid line 
represents the hydrodynamic prediction l/Vl IB3] . The mul- 
tipolarity of the destabilizing mode, see Eq. (nK7 is depicted 
in the inset as a function of temperature. The critical tem- 
perature r° « 45 nK for the parameter values used. 



FIG. 2: Densities of the condensate, the noncondensate, and 

the destabilizing surface mode at T = 40 nK( ) and 10 nK 

( — ) as functions of the radial distance from the symmetry axis 
of the trap. Only part of the condensate boundary is seen in 
the scope of the figure. For clarity, the densities computed at 
T = 10 nK are scaled up by a factor of 10. 



FIG. 3: Thermodynamic (•) and minimum surfon-related (o) 
critical frequencies for vortex nucleation as functions of tem- 
perature. If the rotation of the noncondensate is neglected 
(■), the thermodynamic stability of the vortex state increases 
as the critical temperature is approached from below. Also 

shown is the critical condensation temperature ( ) and the 

frequency for local stabilization of the anomalous vortex core 
mode obtained within the Bogoliubov approximation (BA) 
(— )• 
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